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ABSTRACT 

The spatial clustering amplitude (B gq ) is determined for a sample of 44 powerful 
AGN at z ~ 0.2. No significant difference is detected in the richness of the cluster 
environments of the radio-loud and radio-quiet sub-samples, both of which typically 
inhabit environments as rich as Abell Class ~ 0. Comparison with radio luminosity- 
matched samples from Hill & Lilly (1991) and Wold et al. (2000a) suggests that 
there is no epoch-dependent change in environment richness out to at least z > 0.5 for 
either radio galaxies or radio quasars. Comparison with the APM cluster survey shows 
that, contrary to current folklore, powerful AGN do not avoid rich clusters, but rather 
display a spread in cluster environment which is perfectly consistent with being drawn 
at random from the massive elliptical population. Finally, we argue that virtually all 
Abell class ~ clusters contained an active galaxy during the epoch of peak quasar 
activity at z ~ 2.5. 
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1 INTRODUCTION 



The cluster environments of active galactic nuclei (AGN) 
provide important information for improving our under- 
standing of several aspects of the AGN phenomenon. Over 
the last fifteen years a considerable amount of effort has 
been invested in studying the environments of several dif- 
ferent types of AGN, including : radio-quiet and radio-loud 
quasars (eg. Yee & Green 1984, 1987, Smith, Boyle & Mad- 
dox 1995, 2000, Hall & Green 1998), radio galaxies (Prestage 
& Peacock 1988, 1989, Hill & Lilly 1991), BLLacs (Wurtz et 
al. 1997) and Seyfert galaxies (de Robertis & Yee 1998). In 
addition to its potential for shedding light on the processes 
by which dormant black holes are triggered into AGN, the 
study of clustering environments is an invaluable test of the 
viability of proposed unification schemes (eg. Wurtz et al. 
1997) and models of quasar evolution (eg. Ellingson, Green 
& Yee 1991). In this paper we use hst images to study the 
immediate environments of AGN drawn from well-matched 
samples of powerful radio galaxies (RG), radio- loud quasars 
(RLQ) and radio-quiet quasars (RQQ) at z = 0.2 to explore 
what constraints can be placed on the origin of radio loud- 
ness, the viability of radio loud unification, and the physical 
origin of the dramatic evolution of the quasar population 
between z = 2.0 and the present day. 

Several authors have reported a difference in the rich- 
ness of clustering around radio-loud and radio-quiet quasars 
(eg. Yee & Green 1984, 1987, Ellingson, Green & Yee 1991). 



These studies have shown that, with substantial overlap, 
low to moderate redshift radio-loud quasars are found in 
Abell 0/1 clusters, while radio-quiet quasars rarely inhabit 
clusters as rich as Abell 0. This apparent change in quasar 
environment with radio-power was at least consistent with 
the traditional picture in which RLQs have elliptical host 
galaxies, and are therefore preferentially found in clusters, 
while RQQs inhabit the poorer environments associated 
with lower luminosity disc-dominated Seyfert galaxies. 

However, the results of our own HST host galaxy imag- 
ing programme (McLure et al. 1999, Dunlop et al. 2000), to- 
gether with other recent HST and ground-based studies (eg. 
Hooper et al. 1997, Boyce et al. 1998, Bahcall et al. 1997, 
Schade et al. 2000) show that, at the very least, a significant 
fraction of the RQQ population are located in elliptical host 
galaxies. In fact, the results of our own study lead to the 
conclusion that all objects which could be classified as true 
quasars (i.e. My < —2 3) will have bulge-dominated host 



galaxies. In Section 4.1 we therefore examine the question 
of whether our target RQQs and RLQs, drawn from sam- 
ples well matched in terms of AGN luminosity, are actually 
occupying significantly different cluster environments. 

The other question investigated in this paper is that of 
the evolution of the cluster environments of powerful, radio- 
loud AGN. With respect to radio galaxies, it was shown by 
Prestage & Peacock (1988, 1989) that at low z the envi- 
ronments of powerful frii sources are poorer than their fri 
counterparts, although the substantial overlap between the 
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two populations led Prestage & Peacock to conclude that 
environment did not solely determine radio luminosity. In 
contrast, at higher redshifts, z ~ 0.5, both Yates, Miller 
& Peacock (1989) and Hill & Lilly (1991) found that ra- 
dio galaxies inhabited richer clusters than at low redshift, 
leading Hill & Lilly to conclude that the cluster environ- 
ments of radio galaxies were epoch-dependent. Similar con- 
clusions were reached by Ellingson, Green & Yee (1991) from 
their study of the environments of RLQs. Ellingson et al. 
found that the environments of RLQs also become richer at 
2 > 0.6, although the existence of this epoch-dependence 
has recently been questioned by Wold et al. (2000a). 

In contrast to what might be expected from the results 
outlined above, the raw hst images from our sample of RGs 
and RLQs (Dunlop et al. 2000) show many of the objects to 
be apparently residing in clusters of moderate richness, an 
impression which is confirmed by the a naly sis presented in 
Section^. Motivated by this, in Sections |5.l| fc ^Jwe present 
a comparison between the cluster results for our RG and 
RLQ samples with those for optical and radio luminosity- 
matched sub-samples of the objects studied by Hill & Lilly 
(1991) and Wold et al. (2000a), in order to critically re- 
examine the evidence that RGs and RLQs inhabit environ- 
ments which are a function of cosmic epoch. In Section 6 
we investigate whether there is any evidence in our data for 
a link between AGN properties and the large-scale environ- 
ments of their hosts. Finally, in Section |^ we compare the 
distribution of cluster richness determined for our sample 
with that of the APM cluster survey (Dalton et al. 1997, 
Croft et al. 1997) in light of the constraint imposed by our 
host galaxy study that only galaxies with Mr < —23 are 
capable of producing powerful AGN. Unless otherwise spec- 
ified all cosmological calculations performed in this paper 
assume a cosmology of Ho = 50, qo = 0.5 and A = 0. 



al. (1997) which were not included in our host galaxy study. 
Although image saturation and emission line contamination 
make the Bahcall et al. data flawed with respect to detailed 
host galaxy analysis, these problems are not a concern for 
an environmental study. The data for these objects (3 RLQs 
& 9 RQQs) is of practically identical depth to our own and, 
with the exception of 3C273, falls in the same region of the 
Mv — z plane. 

2.1 Observations and data reduction 

The observations for all 44 objects considered in this paper 
were taken with the HST Wide- Field and Planetary Camera 2 
(WFPC2) . The observations for the 33 objects from our z = 
0.2 host galaxy study utilised the F675W filter, which closely 
approximates the standard Cousins i?-band (Holtzman et al. 
1995), while the data for the objects from the Bahcall et al. 
sample were taken with the F606W (Wide V) filter. Both 
data sets were imaged on the wide-field chips of WFPC2 
(WF2 for the F675W data and WF3 for the F606W data) 
and have a plate-scale of 0.l"/pix. 

The basic data reduction (bias removal & flat-fielding) 
was performed by the HST pipeline. Subsequently, the indi- 
vidual exposures of each source (3 x 600s for the F675W 
data and typically 1100s+600s for the F606W data) were 
then combined and cleaned of cosmic rays using standard 
iraf tasks. The signal-to-noise levels in the final deep im- 
ages correspond to a Cousins f?-band la sensitivity limit of 
fin ~ 26.7 mag arcsec -2 , where the conversion from F606W 
to Cousins i?-band for the Bahcall et al. data assumes a 
typical colour for a z — 0.2 E/Sab galaxy of R c — 606 — 0.3 
(Fukugita et al. 1995). 



3 DETERMINING CLUSTER RICHNESS 



2 THE SAMPLE 

The low-redshift sample studied in this paper comprises 44 
objects and is a combination of that investigated in our 
HST study of z ~ 0.2 AGN host galaxies (McLure et al. 
1999, Dunlop et al. 2000), with data of similar quality from 
the HST archive. The sample investigated in our host-galaxy 
programme consisted of 33 objects, divided into three sub- 
samples of 10 radio galaxies (RG), 10 radio-loud quasars 
(RLQ) and 13 radio-quiet quasars (RQQ). The unique fea- 
ture of this sample is that the two quasar sub-samples were 
originally chosen to be matched in terms of their distribution 
in the optical luminosity-redshift plane (Mv — z) (Dunlop 
et al. 1993) , with the RG and RLQ sub-samples matched in 
both the radio power-redshift plane (Psghz — z), and the 
radio power-spectral index plane (Psghz ~ a) (Taylor et 
al. 1996). Although the selection criteria used in choosing 
these samples were designed to determine the role played 
by host galaxies in the radio-loudness dichotomy and radio- 
loud unification, they are equally valid for investigating the 
nature of their respective cluster environments. One of the 
objects from this sample, the RQQ 1549+203, has however 
been excluded from this analysis due to the known presence 
of a foreground cluster (Stocke et al. 1983). In order to in- 
crease the numbers of objects studied we have also included 
in our sample 12 quasars from the HST study of Bahcall et 



The method adopted in this study for quantifying the rich- 
ness of the AGN environments was to determine their re- 
spective spatial clustering amplitudes (B gq ) (Longair & 
Seldner 1979). This is a standard technique which, although 
originally designed to investigate clustering around radio 
galaxies, has subsequently been successfully applied to the 
environments of a wide variety of active and inactive galax- 
ies (eg. Yee & Lopez-Cruz 1999, Ellingson, Yee & Green 
1991, Prestage & Peacock 1988). A detailed description of 
the derivation of B gq is given in Longair & Seldner (1979) 
and consequently only a brief outline is provided here. The 
first stage in the calculation is the determination of the an- 
gular correlation function, defined as : 



n (e) sn = N g [1 + w (9)] sn 
where 



(1) 



w(9)=A gq 9 1 -' (2) 

and A gq quantifies the excess in the number of galaxies in 
the vicinity of the source as compared to the expected back- 
ground contribution N 9 . Provided that < 1 the value of 
A gq can be directly calculated from the data using the ex- 
pression: 



Agq — 



^ - 1 



(3) 
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where N t is the total number of galaxies counted within a 
radius of 9 radians from the target (excluding the target 
itself), and Nb is the expected number of background galax- 
ies within the same radius. In order to directly compare the 
clustering around objects which cover a range of redshifts it 
is then necessary to de-project this angular correlation into 
its spatial equivalent which is defined by: 



n (r) SV = p a [1 + e (r)] SV 
where 

e{r) = B aq r-^ 



(4) 



(•») 



and B gq is the desired spatial clustering amplitude. By in- 
voking the simplifying assumption that the clustering of 
galaxies is spherically symmetric around the central object it 
can be shown (Longair & Seldner 1979) that the spatial and 
angular clustering amplitudes are related by the expression: 

-3 



So 



AggN g 
Ij(j> (z) 



D 



1 + z 



(6) 



where D is the effective angular diameter distance to the tar- 
get and 4> ( z ) is the integrated field-galaxy luminosity func- 
tion at the redshift of the target. The quantity L t is a con- 
stant which has a value of 3.78 for the canonical field-galaxy 
value of 7 = 1.77 (Groth & Peebles 1977). 

For a galaxy appearing on the same WF CCD as one of 
the AGN in this sample to be counted as a possible cluster 
member it had to satisfy two criteria. Firstly, its projected 
distance from the central object had to be less than the 
counting radius; defined as the distance of the central ob- 
ject from the nearest edge of the CCD. At the sample median 
redshift of z=0.2 this radius corresponds to a projected met- 
ric radius of typically 180 kpc. Although this is undoubtedly 
a small counting radius compared to the usual 500 kpc or 
1 Mpc counting radii adopted in most cluster studies, pro- 
vided that the clustering of galaxies around the AGN does 
have a slope of 7 = 1.77, the restriction of only having in- 
formation of the central regions of the clusters should not 
prevent a reliable determination of the enhancement of as- 
sociated galaxies relative to the field population. Indeed, as 
will be discussed in Section 3.2, there is good evidence that 
the clustering around the AGN studied here does follow the 
expected slope of 7 = 1.77, at least on average. Further- 
more, in their recent study of low- z Abell clusters, Yee & 
Lopez-Cruz (1999) found that reducing the counting radius 
from 1 Mpc to 200 kpc only altered the spatial clustering 
amplitude determination at the 10 percent level. 

The second selection criterion was that galaxies had to 
lie in the magnitude interval m, — 1 — > + 2, where m* is 
the apparent magnitude corresponding to M* at the AGN 
redshift. This magnitude interval is a compromise which is 
designed to include those cluster galaxies which contain the 
majority of the cluster mass. The faint limit of m* + 2 should 
probe deep enough into the cluster luminosity function to 
be sensitive to any enhancements in galaxy density, without 
risking the possibility of missing clusters due to the back- 
ground galaxy counts rising more rapidly than the faint-end 
of the cluster luminosity function at m > m* + 2. The bright 
limit of nu — 1 is set to avoid problems associated with small 
number statistics, although relaxation of this limit to m* — 3 
has a negligible effect on the results presented in Section ^ 
An additional advantage with this choice of magnitude in- 
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-22.32 
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0.50 


-> 0.75 


0.0078 


-22.11 


0.50 
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-> 2.00 


0.0068 


-22.52 


1.28 



Table 1. The Schechter luminosity function parameters used to 
predict the background galaxy counts. Column 1 lists the four 
redshift bins adopted. The parameters for the lowest redshift bin 
have been taken from Yee & Lopez-Cruz (1999). The parameters 
in the other three redshift bins have been taken from Lilly et al. 
(1995) with a conversion of M* deduced assuming typical colours 
of B - B AB = 0.17 (Metcalfe et al. 1991) and R - B = 1.45 
(Fukugita et al. 1995). 
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Figure 1. A schematic of the Wide Field and Planetary Camera 
2 showing the spatial coverage at the median redshift of z = 0.2. 
The 32 objects imaged during the AGN host galaxy study of 
McLure et al. (1999) and Dunlop et al. (2000) are centred on 
WF2. The 12 objects included from the sample of Bahcall et 
al. (1997) were imaged on WF3. The extra separation of WF2 
and WF4 allowed the distribution of background galaxies and 
the slope of the correlation function to be investigated (Section 

I- 



terval is that, due to the host galaxies of the AGN studied 
in our hst imaging programme having average luminosities 
one magnitude brighter than M* (Dunlop et al. 2000), this 
interval typically corresponds to the m s -^m 9 + 3 range 
adopted for the Abell-type calculation of cluster richness 
around z ~ 0.5 radio galaxies by H ill & Lilly (1991). The 
comparison presented in Section 5.1 of the clustering around 
the z ~ 0.2 radio galaxies studied here, with the results ob- 
tained by Hill & Lilly (1991), can consequently be performed 
in a more transparent manner. 



3.1 Galaxy counts and the luminosity function 

A crucial element in the calculation of the spatial cluster- 
ing amplitude is the form of the luminosity function used in 
the normalization of equation ^. Given that the derivation of 
equation ^ is dependent on the assumption that the observed 
background galaxy counts can be predicted from integrat- 
ing the galaxy luminosity function along the line of sight, 
it is essential that the luminosity function chosen should 
at least be consistent with the background counts. In or- 
der to ensure this we have adopted a similar approach to 
that employed by Wold et al. (2000a), by considering the 
form of the luminosity function in four redshift bins z=(0.0, 



© 0000 RAS, MNRAS 000, 000-000 



4 R.J. McLure 



0.2), (0.2, 0.5), (0.5, 0.75) & (0.75, 2.0). The Schechter func- 
tion parameters of the four luminosity functions are listed in 
Table [j]. In the lowest redshift bin we have adopted the lu- 
minosity function determined by Yee & Lopez-Cruz (1999) 
from their J?-band study of Abell clusters in the redshift 
range 0.02 < z < 0.18. The characteristic magnitude and 
faint-end slope of this Schechter function are nearly identi- 
cal to those determined for the field-galaxy population by 
Loveday et al. (1992), but with a normalization some 30% 
higher. Due to the fact that the Loveday et al. determina- 
tion from the Stromlo-APM survey is dominated by objects 
with z < 0.1, it was felt that the Yee & Lopez-Cruz lumi- 
nosity function was a fairer representation of the z ~ 0.2 
galaxy population. The adoption of the higher normaliza- 
tion also has the advantage that we are correspondingly less 
likely to be overestimating the spatial clustering amplitudes. 
In the three highest redshift bins the parameters are taken 
from Lilly et al. (1995), with the appropriate conversion of 
the characteristic magnitudes from the Mab system to the 
Cousins i?-band (see Table 0). The predicted number counts 
in a particular magnitude range (mi,rri2) were then calcu- 
lated by integrating the following function: 



Blank Fields 



AGN Fields 



(m, z)5m 



/SV 
{ 8z 



8z 



(7) 



where the limiting redshift of z — 2 was chosen since, in 
the magnitude range investigated here (R < 24), the con- 
tribution from background galaxies at z > 2 is negligible. 
A comparison between the predicted and measured galaxy 
counts for both the target (WF2+WF3) and control CCDs 
(WF4) is shown in Fig |2j. It can be seen that the predicted 
background galaxy counts (left panel) agree well with those 
measured in the magnitude range (21 < R < 23), where 
the control chips should be relatively unaffected by contam- 
ination from cluster members. However, in the magnitude 
range 18.5 < R < 20.5 there is an excess of detected back- 
ground galaxies over that predicted, as expected given that 
at our median redshift of z = 0.2 this range corresponds to 
M* — > M* + 2. The predicted background counts are in good 
agreement with those determined by both Yee & Lopez-Cruz 
(1999) and Metcalfe et al. (1991) and were subsequently used 
to remove the background contribution in the calculations 
of the spatial clustering amplitude. The comparison between 
the predicted background counts and the counts from the 
AGN target CCDs shows qualitatively that inside a radius 
of ~ 200 kpc we are clearly detecting a large excess of cluster 
galaxies. 



3.2 The slope of the correlation function 

The usual procedure when determining the richness of a 
cluster environment via its spatial clustering amplitude is 
to assume that the slope of the two-point correlation func- 
tion has the canonical value of 7 = 1.77 (Groth & Peebles 
1977). Here we investigate whether this is in fact a reason- 
able assumption for the objects in this sample. As a result 
of the relative positions of the WF2 and WF4 CCDs on the 
WFPC2 (see Fig [l]) it is possible to perform galaxy ring 
counts out to a projected radius of ~ 800 kpc for the 32 ob- 
jects imaged through the F675W filter, under the assump- 
tion of spherical symmetry. Fig ^| shows the average radial 
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Figure 2. The left-hand panel shows the background galaxy 
number counts from the 32 WF4 observations from McLure et 
al. (1999) and Dunlop et al. (2000), at an average separation 
from the target AGN of ~ 500 kpc. The solid line is the predicted 
number counts from integrating the galaxy luminosity function 
along the line of sight (see text for discussion). The right-hand 
panel shows the galaxy counts from the target WF chips of all 44 
objects at a separation of < 200 kpc. A clear excess of galaxies is 
obvious at R < 23. 
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Figure 3. Radial profile showing the average ratio of excess clus- 
ter galaxies (N c ) to background galaxies (N g ). The width of the 
bins changes from 60 kpc to 120 kpc at a radius of 220 kpc in order 
to maintain comparable signal-to-noise. Also shown is the best-fit 
power-law correlation function (dotted line). The best fit value of 
7 is l.SSj^S'Jg where the errors refer to the A\ 2 — ^X^in = ^ 
confidence limit. 



profile of the ratio of cluster galaxies to background galax- 
ies (N c /Ng) measured in metric aperture bins. In order to 
keep the signal-to-noise level of the bins approximately con- 
stant their width is increased from 60 kpc to 120 kpc at a 
radius of 200 kpc. The dotted line in Fig | shows the min- 
imum x 2 fit f° r the angular correlation function, as defined 



in equation 0, which has a value of 7 = 1.. 



It appears 



therefore that, at least on average, the form of the correla- 
tion function for the objects in this sample is consistent with 
the standard value of 7 = 1.77. Consequently, we adopt this 
standard value in our calculation of the spatial clustering 
amplitudes in order to allow easy comparison with results 
in the literature. 
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Figure 4. Histograms showing the distribution of clustering am- 
plitudes displayed by the three sub-samples. Also shown in the 
approximate level of the clustering amplitudes of field galaxies 
(dotted line), here taken as B 9g = 60 



Sample 


N 


< B gq > 


Median 


All 


44 


365 ± 62 


241 


RG 


10 


575 ± 165 


321 


RLQ 


13 


267 ± 51 


247 


RQQ 


21 


326 ± 94 


209 



Table 2. Summary of the spatial clustering amplitude results. 
Column two details the number of objects in each sub-sample. 



4.1 Do RLQs reside in richer environments than 
RQQs? 

The clustering results for the two quasar sub-samples shown 
in Tables 2 & ^ indicate that there is no significant dif- 
ference in the richness of the cluster environments of the 
RQQs and RLQs studied here. It can be seen from Table 
2 that the two sub-samples display mean and median clus- 
tering amplitudes which are statistically consistent, with the 
Kolmogorov-Smirnov (KS) test returning a significance level 
of only p=0.54. 

The finding that RLQs inhabit environments with spa- 
tial clustering amplitudes of the order ~ 200 — > 300 Mpc 1,77 
is in good agreement with previous studies. For example, 
the large-scale study of Ellingson, Yee & Green (1991) found 
that at z < 0.4 the average clustering amplitude of RLQs 
was 210 ± 70 Mpc 1,77 . Moreover, the clustering study per- 
formed by Fisher et al. (1996) using the hst data for the 
20 quasars in the Bahcall et al. (1997) sample (all of which 
are common to the 44 object sample studied in this paper) 
found the 6 RLQs in that sample to have an average clus- 
tering amplitude of ~ 256 ± 55 Mpc 1 ' 77 . 

In contrast, the finding that the environments of RQQs 
are indistinguishable from those of RLQs is not in general 
agreement with the literature, with the majority of previ- 
ous studies concluding that RQQs tend to inhabit poorer 
environments. For example, the studies of both Ellingson, 
Yee & Green (1991) and Smith, Boyle & Maddox (1995) 
find the clustering around RQQs at z < 0.3 to be perfectly 
consistent with that of field galaxies, B gq /B gg = 1.0, while 
in contrast even the median value found here corresponds 
to B gq / B gg = 3.8, taking B gg — 60. The exact cause of this 
discrepancy is difficult to determine due to the large number 
of factors which influence the calculation of B gq . One possi- 
ble explanation for the disagreement between these results 
and those of Smith et al. is that ~ 80% of their objects are 
fainter than My = —23, which could explain the detection 
of environments typical of Seyfert galaxies (eg. de Robertis 
& Yee 1998). This issue was looked at by Smith et al., who 
found no significant difference in their results if they divided 
their objects into high and low luminosity sub-samples, al- 
though this split would still have resulted in ~ 50% of the 
objects in the high-luminosity sub-sample being fainter than 
M v = -23. 

It is worth noting at this point that the results pre- 
sented by Fisher et al. (1996) of the clustering around the 
quasars in the Bahcall et al. (1997) sample do agree with the 
results presented here. Fisher et al. found an average value 



4 RESULTS 

The results of the B gq calculations for the three sub-samples 
are listed in column 6 of Table ^| and are shown graphically 
in Figs ^ & ^ In this section we discuss the results for the 
separate sub-samples in the context of the results of our 
host galaxy study (McLure et al. 1999, Dunlop et al. 2000) 
and previous stu dies of similar objects from the literature. 
In Sections 5.1 & 5.2 we proceed to perform a more detailed 



comparison between our results and those obtained by Hill 
& Lilly (1991) and Wold et al. (2000a) at higher redshift. 



data for all of which are included in the 21 RQQs stud- 
ied here. Given that both of our studies are based on HST 
data this obviously raises the concern that the small field 
of view of WFPC2 has led to an overestimate of the rich- 
ness of the RQQs environments. However, there are several 
reasons for believing that this is not a serious issue. Firstly, 
as demonstrated by Fig ^, there is no real suggestion that 
the enhancement in associated galaxies is falling off faster 
than 7 = 1.77. Secondly, given the good agreement between 
our results for the RLQ sub-sample and previous studies, 
the form of clustering around the RQQs would have to be 
substantially different from that around RLQs in order to 
account for the difference. Although this is undoubtedly a 
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possibility, it is not supported by our data, with the F675W 
images showing no tendency for the RQQs to have fewer as- 
sociated galaxies at radii of 200 < r < 700 kpc. It appears 
therefore that the similarity in environments of our RQQ 
and RLQ samples probably results from the close match- 
ing of the samples in the Mv — z plane. This conclusion is 
supported by the preliminary results of Wold et al. (2000b), 
who also find no difference in cluster environment between 
their optically matched samples of RQQs and RLQs in the 
redshift range 0.5 < z < 0.8. 



4.2 Are the environments of RLQs and RGs 
consistent with unification? 

With the results of our host-galaxy study showing that our 
RG and RLQ samples are virtually identical in terms of 
scalelength, luminosity, and R — K colour, it was expected 
that their cluster environments would also be indistinguish- 
able. The results presented in Table ^ and Fig ^ broadly 
confirm this expectation, although it is clear that the clus- 
ter environments of the two samples are not as similar as 
suggested by the optical properties of their host galaxies. 
The application of the KS test shows the two distributions 
to be distinguishable only at the la level (p = 0.25), a differ- 
ence which is clearly due to the three RGs with B gq > 1000 
Mpc 1 ' 77 . Considering the large errors associated with the 
B gq calculation, and the small number statistics, perhaps 
the strongest statement that can be made is that these re- 
sults present no problem to the orientation-based unification 
of RGs and RLQs due to differences in cluster environment. 
The possibility of a correlation between environment and 
radio power is investigated in Section H 



4.3 Abell classification 

In this section we attempt to provide a transformation be- 
tween the spatial clustering amplitudes and the traditional 
Abell cluster classification (Abell 1958, Abell, Corwin & 
Olowin 1989). Unfortunately, due to numerous sources of 
systematic error associated with Abell counts, the correla- 
tion between Abell class and B gq has a large scatter, and 
a consensus on the appropriate calibration has not been 
reached in the literature. In their study of 47 z < 0.2 
Abell clusters Yee & Lopez-Cruz (1999) proposed a lin- 
ear scheme whereby adjacent Abell classes are separated 
by AB gq = 400 Mpc 1 ' 77 , with a normalization such that 
Abell class clusters have a spatial clustering amplitude of 
Bgq = 600 Mpc 1 ' 77 . Here we have chosen to classify our clus- 
ter measurements using a scheme which is identical to that 
proposed by Yee & Lopez-Cruz save for a re-calibration such 
that Abell class clusters correspond to B gq — 300 Mpc 1,77 
(Table g). 

The main reason for this re-calibration is that, unlike 
the clusters of class > 1, the Abell clusters studied by Yee 
& Lopez-Cruz were not selected randomly, but specifically 
chosen as being rich. Consequently, it is almost certainly the 
case that a figure of B gq — 600 Mpc 1,77 is not representative 
of Abell clusters as a whole. Moreover, spatial clustering 
amplitudes of B gq = 300 & 700 Mpc 1 ' 77 for Abell classes 
& 1 are much more representative of the findings of sev- 
eral previous studies (eg. Longair & Seldner 1979, Prestage 
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Table 3. Sample details and spatial clustering results. Column 
5 details the 5GHz radio luminosity (where known) in units of 
WHz" 1 sr _1 . The errors have been calculated using the con- 
servative prescription of Yee Sz Lopez-Cruz (1999): -g- 2 ^ = 

[(N t -N b )+1.3 2 N b } 1/2 
N t -Nb 



k. Peacock 1988, Andersen & Owen 1994). Evidence that 
the calibration proposed here is also reasonable for the very 
richest clusters comes from the fact that the four richest 
clusters studied by Yee & Lopez-Cruz had a mean value of 
Bgq — 2225 Mpc 1 ' 77 , here corresponding to Abell class 4/5, 
with no clusters found with B gq > 2300 Mpc 1,77 It can be 
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Figure 5. The 3 gq — z distribution for the three sub-samples. Also shown are the approximate Abell cluster classifications according to 
the linear scheme shown in Table 4. 
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Table 4. The proposed relation between spatial clustering am- 
plitude and Abell classification. This scheme is identical to that 
proposed by Yee & Lopez-Cruz (1999) but with Abell Class 
clusters normalized to B gq = 300 instead of 700 Mpc 1,77 (see 
text) 



seen that using our calibration of the B gq — Abell transfor- 
mation the results presented in this section imply that the 
average clustering around our AGN corresponds to Abell 
class ~ 0, with the large scatter ranging from several objects 
consistent with no galaxy enhancement, to several clusters 
of Abell class 2/3. 



5 DO THE ENVIRONMENTS OF 

RADIO-LOUD AGN CHANGE AT Z=0.5? 

The clustering results for the radio galaxy and radio-loud 
quasars sub-samples presented in the previous section offer 
an opportunity to re-examine the evidence for an epoch- 
dependent change in cluster richness at z ~ 0.5. In this 
section we investigate this issue by comparing our results 
with those published for radio galaxies by Hill & Lilly (1991) 
and the recent work on radio-loud quasars by Wold et al. 
(2000a). 

5.1 Radio galaxies 

In their study of the environments of radio galaxies Hill & 
Lilly (1991) studied four groups of ~ 10 objects at z ~ 0.5, 
each of which were equally spaced in radio power from Log 
P2GHz = 23 — + 27 WHz -1 . To compare their results with 
those for our radio galaxy sample we selected objects from 
their two medium radio luminosity groups, subject to the 
redshift constraint 0.40 < z < 0.53, which was chosen to 
match the range in redshift displayed by our radio galaxies 
(0.12 < z < 0.26). This selection procedure produced 16 ob- 
jects, two of which were quasars and subsequently rejected 
to leave a final sample of 14 objects. The P2GHZ— Ms and 
P2GHZ — z distributions of the Hill & Lilly sub-sample and 



our radio galaxy sample are shown in the left and middle 
panels of Fig ^| The two samples are clearly well matched, 
with the two-dimensional KS test (2DKS) showing their re- 
spective P2&HZ— Ms distributions to be statistically indis- 
tinguishable (p = 0.2). 

To facilitate a comparison between the cluster richness 
determinations of Hill & Lilly and our own results it was 
necessary to transform between B gq and their A*o.5 values. 
The value of A'o.s is the background subtracted number of 
galaxies within a radius of 0.5 Mpc, in the magnitude range 
m g — » m g + 3, where m g is the apparent magnitude of the 
radio galaxy. As previously mentioned in Section this 
magnitude range should be well matched to our choice of 
nu — 1 -^m* + 2, and to make the transformation we have 
used Hill & Lilly's own determination of B gq — 30No.$. The 
resulting B gq — z distributions for the two samples are shown 
in the right-hand panel of Fig ^. It can be seen that the two 
distributions are very similar, an impression confirmed by a 
KS test (p = 0.77), with no indication that the members of 
the low- 2 sample inhabit poorer environments, or that they 
display a smaller range in environment richness. Although 
we are dealing with small samples, this result demonstrates 
that, at least for these two well matched samples, there ap- 
pears little evidence for a epoch-dependent change in FRII 
radio galaxy environments at z ~ 0.5. In light of this, we 
now move on to look for a change in the environments of 
radio- loud quasars with redshift. 

5.2 Radio-loud quasars 

The recently published study of 21 RLQs in the redshift 
interval 0.5 < z < 0.8 by Wold et al. (2000a) provides a good 
opportunity for comparison with the results from our z ~ 
0.2 RLQ sample. Wold et al. also classify the environments 
of their quasars via the spatial clustering amplitude and 
their results can therefore be directly compared to our own. 
Taken as a whole the results for their sample are in excellent 
agreement with our own, with their mean figure of B gq = 
265±65 Mpc 1,77 being almost identical to our mean of B gq = 
267 ± 51 Mpc 177 . 

To disentangle the effects of redshift and radio lumi- 
nosity we have excluded the 7 most radio luminous of the 
Wold et al. sample (Log P408 > 26.5 WHz _1 sr _1 ) to leave a 
sub-sample of 14 objects which are well matched to our 12- 
object RLQ sample. The distribution of the two samples on 
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Figure 6. The left-hand and middle panels show the matching of our RG sample (open circles) with the sub-sample of the objects 
studied by Hill & Lilly (1991) discussed in the text (filled circles). The right-hand panel shows a comparison between their respective 
spatial clustering amplitudes. 
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Figure 7. The left-hand and middle panels show the matching between our RLQ sample (open circles) and the sub-sample of the objects 
studied by Wold et al. (2000a) discussed in the text (filled circles). The right-hand panel shows a comparison of their respective spatial 
clustering amplitudes. 



the P408— Mb and P408 — z planes is shown in the left-hand 
and middle panels of Fig The two samples can be seen to 
be well matched, with the 2DKS test returning a probability 
of p — 0.51 that the two samples are drawn from the same 
P408— Ms distribution. From the comparison of the spatial 
clustering amplitudes of the two samples shown in the right- 
hand panel of Fig (?|, there is again no suggestion that the 
high-z sample displays systematically higher values of B gq , 
although it is noticeable that the scatter in B gq is somewhat 
higher in the high-z sample. This is in fact the conclusion 
arrived at by Wold et al., who also found no evidence for an 
epoch-dependent change in RLQ environments when com- 
paring their results with those of Ellingson, Yee & Green 
(1991) at low redshift. 



6 EVIDENCE FOR A LINK BETWEEN AGN 
PROPERTIES AND HOST CLUSTER 
ENVIRONMENT 

In this section we explore whether there is any evidence in 
our data for a link between AGN properties and the large- 
scale environments of their hosts. This issue is of interest 
because it offers the potential to determine whether the ra- 
dio and optical luminosity of powerful AGN are primarily 
influenced by large-scale environmental factors such as den- 



sity of the ICM and galaxy interactions, or alternatively, 
whether the physical properties of the active nucleus itself, 
such as black-hole mass and accretion rate, are the dominant 
influence. 

The three panels shown in Fig |^ display the relationship 
between B gq and respectively, estimated central black-hole 
mass, quasar accretion efficiency and radio luminosity. All 
three panels feature only those objects from the sample im- 
aged in the McLure et al. (1999), and Dunlop et al. (2000) 
host-galaxy study in order to produce a homogeneous data- 
set. The objects from the Bahcall et al. (1997) sample have 
been excluded on the basis that, due to the saturated nature 
of their data, reliable nuclear luminosities are unavailable. 

In the left-hand panel of Fig |^ we have plotted B gq 
against black-hole mass, as calculated by a straight applica- 
tion of the galaxy bulge mass-blackhole mass relation pub- 
lished by Magorrian et al. (1998) to the host galaxy mod- 
elling results presented in Dunlop et al. (2000). A general 
trend for black-hole mass to increase with spatial clustering 
amplitude can be seen, although the correlation is not sta- 
tistically significant, with a Spearman rank test returning a 
probability of p — 0.38 of no correlation. The clearest sig- 
nal comes from the ten radio galaxies, although even here 
Kendall's tau test still returns only a marginally significant 
result; p — 0.09. However, the lack of a clear correlation 
between B gq and black-hole mass is perhaps not entirely 
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Figure 8. Plots of the B gq versus black-hole mass (left), quasar nuclear luminosity as a fraction of the Eddington limit (middle) and 
radio luminosity at 408MHz (right). In all three plots RQQs are represented by filled circles, RLQs are represented by open circles and 
RGs are represented by open squares. 



unexpected. Even if the basic form of the cluster luminosity 
function does not change significantly with cluster richness, 
there will be a corresponding increase in the number of lumi- 
nous, bulge-dominated galaxies which are capable of hosting 
powerful AGN as cluster richness increases. Consequently, 
in rich clusters, it is not only the central brightest-cluster 
galaxy which is a potential AGN host, but also the next few 
lower-ranked cluster galaxies. This then will inevitably in- 
troduce significant horizontal scatter into the B 9 q-blackhole 
mass relation, regardless of the 0.5 dex scatter in the Magor- 
rian relation itself. 

The middle panel of Fig ^ shows the relationship be- 
tween B gq and the nuclear luminosity of the quasars, ex- 
pressed as a fraction of their Eddington limit, under the 
assumption of a central black-hole mass as predicted by the 
Magorrian relation. As with the previous plot, a trend can 
be seen for the Eddington fraction to increase with cluster 
richness, although again this is not statistically significant, 
with a rank test returning a probability of p = 0.5 that no 
correlation is present. However, it is perhaps interesting to 
note that the three RQQs which are radiating closest to their 
predicted Eddington limit are also ranked first, second and 
fourth in terms of their cluster richness. 

Finally, the right-hand panel of Fig ^ examines whether 
there is a correlation between the radio luminosity of the 
radio-loud AGN in our sample, and the density of the en- 
viroment into which their radio jets are expanding. As with 
the other two plots in Fig |^ a positive correlation is again 
suggested to the eye, but is once more found not to be signif- 
icant, with a rank test returning a probability of p = 0.47 of 
no correlation, even after the removal of the obvious outlying 
radio galaxy. It is also noteworthy that the apparent positive 
trend between radio luminosity and B gq is entirely defined 
by the ten radio galaxies, with the RLQs showing no evi- 
dence for a connection between radio luminosity and cluster 
richness. Consequently, at least for this sample of objects, 
there appears no real evidence that the density of a radio- 
loud AGN's cluster environment plays the dominant role in 
determining its radio luminosity. However, it is also clear 
that the limited dynamic range of this sample, combined 
with the inherently large uncertainties associated with spa- 
tial clustering amplitudes, means that the results presented 
here are still perfectly consistent with a weak correlation 
being present. 




Figure 9. A comparison between our AGN cluster number counts 
with those of the APM cluster survey using a bin size of AB gq = 
400. The solid circles are our uncorrected number counts, while 
the open circles are weighted using the inverse of the mean B gq 
in each bin (see text). The solid line is a power-law with slope - 
3.7 which fits the APM number counts and has arbitrary vertical 
normalization. The vertical dotted line indicates where the APM 
survey becomes incomplete. 



7 COMPARING AGN ENVIRONMENTS TO 
THE APM CLUSTER SURVEY 

In Section |E| it was shown that with samples of RQQs 
and RLQs which are carefully matched in terms of opti- 
cal luminosity there is no evidence that RQQs occupy sys- 
tematically lower-density cluster environments. This result 
therefore makes it of interest to re-examine the question of 
whether quasars actually avoid rich clusters, as sometimes 
claimed, or if alternatively the apparent lack of quasars in 
rich clusters is simply a reflection of the fact that rich clus- 
ters are relatively scarce. 

In order to address this question we have compared our 
cluster results with cluster number counts from the APM 
survey, which are well modelled by a power-law with slope 
~ —3.7 (Dalton et al. 1997). To allow a fair comparison of 
the relative numbers of clusters it is necessary to correct for 
the fact that the number of potential AGN host galaxies in- 
creases with cluster richness. To allow for this effect we make 
the assumption that the galaxy population of all our clus- 
ters can be described by a Schechter function with constant 
characteristic magnitude and faint-end slope, and a normal- 
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ization simply proportional to B gq . Consequently, in Fig g 
we have simply weighted the relative number counts of clus- 
ters by the inverse of the mean B gq in each bin. As can be 
seen in Fig ^ this simple correction naturally produces num- 
ber counts which are perfectly consistent with the power-law 
found by the APM survey. 

A further interesting result arises if one considers the 
space densities of clusters from the APM survey as calcu- 
lated by Croft et al. (1997). With the calibration that in the 
scheme used in the analysis of the APM survey, a cluster 
of richness R=40 corresponds to roughly Abell class (ie. 
B gq ~ 300 Mpc 1,77 ), it is possible to extrapolate the figures 
provided by Croft et al. to estimate that the space density of 
all clusters richer than R> 40 is ~ 10 -5 Mpc -3 . If one allows 
for the fact that up to 50% of quasars could be obscured in 
the optical-UV by a dusty torus, this abundance of clusters 
is in reasonable agreement with the peak space density of 
quasars at z ~ 2.5 (Warren, Hewett & Osmer 1994). 

This line of argument leads to an interesting conclusion 
concerning the fraction of massive cluster galaxies which 
were active at z ~ 2.5. In their photometric study of 83 
Abell clusters Schneider, Gunn & Hoessel (1983) determined 
the absolute magnitudes of the three brightest galaxies in 
each cluster. When converted to our cosmology, the mean 
magnitude of the first ranked galaxy in clusters of Abell 
class corresponds to Mr ~ —23.1 ± 0.5, where a colour 
of r — R — 0.35 has been assumed (Fukugita et al. 1995). 
This is in excellent agreement with the figure of Mr ~ — 23 
which was found to be the necessary bulge luminosity for a 
galaxy to host a powerful (Alv < —23) quasar by Dunlop et 
al. (2000). Furthermore, because the mean luminosity found 
by Schneider, Gunn & Hoessel for the second ranked galaxy 
in Abell class clusters was approximately one magnitude 
fainter than that of the first ranked galaxies, we can conclude 
that the average number of massive elliptical galaxies in a 
Abell class cluster which are capable of hosting a powerful 
quasar is close to unity. Consequently, because of the steep 
decline in the number of clusters with increasing richness 
(Fig ^), the close agreement between the peak space density 
of quasars and that of clusters of Abell class > leads to 
the conclusion that virtually all massive cluster ellipticals 
were active at z ~ 2.5. This conclusion is unaffected by the 
fact that the very richest clusters may have contained 4 or 
5 active galaxies at this epoch because the rarity of such 
clusters means that they have a negligible affect on quasar 
numbers. 



8 CONCLUSIONS 

The cluster environments of a matched sample of 44 pow- 
erful radio-loud and radio-quiet AGN at z ~ 0.2 have been 
analysed using the spatial clustering amplitude method. The 
three main conclusions arising from this study can be sum- 
marized as follows: 

Firstly, we find no evidence to suggest that RQQs in- 
habit poorer cluster environments than RLQs or RGs. All 
three classes of AGN are located in environments which, on 
average, are comparable to Abell Class 0, although there is 
a large amount of scatter. 

Secondly, by comparing our results with those for sam- 
ples from Hill & Lilly (1991) and Wold et al. (2000a) 



which are matched in terms of redshift and radio luminosity, 
we have re-examined the evidence for an epoch-dependent 
change in the environments of RLQs and RGs. We find no 
evidence that the environments of RGs and RLQs become 
significantly richer at z ~ 0.5 as has been reported by Hill & 
Lilly (1991) and Ellingson, Yee & Green (1991) respectively. 

Finally, via comparison with the APM cluster survey we 
conclude that the distribution of AGN cluster environments 
is consistent with having being drawn at random from the 
general cluster distribution. Furthermore, because the clus- 
ter population is dominated by clusters of Abell class ~ 
which, on average, only contain one galaxy capable of host- 
ing a powerful quasar, we argue that the close agreement 
between the space density of clusters and the peak space 
density of quasars at z ~ 2.5 suggests that practically all 
massive cluster ellipticals were active at this epoch. In this 
scenario, the apparent lack of powerful AGN in rich clusters 
at the present day simply reflects the comparative scarcity 
of high density environments. 
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